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Abstract 

Lattice QCD is rapidly progressing toward being able to reliably compute the electric dipole mo- 
ment of the neutron as a function of the strong CP-violating parameter 6. Present day calculations 
are performed at unphysical values of the light quark masses, in volumes that are not exceptionally 
large and at lattice spacings that are not exceptionally small. We use chiral perturbation theory 
to determine the leading contributions to the neutron electric dipole moment at finite volume, and 
in partially-quenched calculations. 
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I. INTRODUCTION 



CP- violation is still a mystery. Current measurements of CP-violating processes in the kaon 
and B-meson sectors would suggest that the single phase in the CKM matrix provides a 
complete description. However, the baryon asymmetry of the universe cannot be described 
by this phase alone, and there are additional sources of CP-violation that await discovery. 
The recent revelation that neutrinos have non-zero masses has presented us with the pos- 
sibility of CP-violation in the lepton sector. With both Dirac and Majorana type masses 
possible, CP-violation in the neutrino sector is likely to be far more intricate than in the 
quark sector. The significant number of experiments operating in, and planned to explore 
the neutrino sector will greatly improve our knowledge in this area in the not so distant fu- 
ture. It has been a puzzle for many years that there is the possibility of strong CP- violation 
arising from the 9 term in the strong interaction sector, but that there is no evidence at this 
point in time for such an interaction. The naive estimate for the size of observables, such as 
the electric dipole moment (edm) of the neutron, induced by such an interaction is orders 
of magnitude larger than current experimental upper bounds, thereby placing a stringent 
upper bound on the coefficient of the interaction, 0. An anthropic argument that compels 
9 to be small does not yet exist and so it is likely that there is an underlying mechanism, 
such as the Peccei-Quinn mechanism and associated axion, that eliminates this operator. 
With the increasingly precise experimental efforts to observe the neutron edm P, , it is 
important to have a rigorous calculation directly from QCD. 

Lattice calculations of the neutron edm j2|, |3|, |5j, |6||, ||7|] in terms of the strong CP-violating 
parameter are continually evolving toward a reliable estimate that can be directly compared 
with experimental limits and possible future observations. The latest generation of lattice 
calculations respect chiral symmetry and lattice spacing effects have been relegated to 0{a^). 
However, the calculations are performed in modest finite volumes and at quark masses that 
are larger than those of nature. In this work we explore the impact of finite volume on such 
calculations and also examine the quark mass dependence of partially-quenched calculations. 

The QCD Lagrange density in the presence of the CP-violating 6'-term is 

C = qilpq - q^m.q^ - q^m\qL - -G^^^^^G^^J + Q^^G^i^G^^-^^^ , (1) 

where G^^^^"' = ie'^^^^Gj^^^ e^^^a ^ ^^^^ ^^lere q = {u,df for two-fiavor QCD. Chiral 
redefinitions of the quark fields modify the coefficient of the GG operator through the strong 
anomaly, and as a consequence it is the quantity 9 = 9 — arg(det(mq)) that has physical 
meaning. For our purposes it is convenient to start with a Lagrange density where niq is 
real and diagonal, and 9 in eq. (0) is equal to 9. One can then remove the GG operator 
by a chiral transformation, qjji — > e^'^^ ^"^ qj r, and qjL — > e~'^'^^f'^qjL subject to the constraint 
that 9 = — 0j- Under this transformation the elements of the quark mass matrix become 

The low-energy effective field theory (EFT) describing the behavior of the pseudo- 
Goldstone bosons associated with the breaking of chiral symmetry is, at leading order. 



C = ^Tr 



D^^ D^St + A ^ Tr [ m,St + mJS J , (2) 
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where / ~ 132 MeV is the pion decay constant, the covariant derivative describing the 
couphng of the pions to the electromagnetic field is -D^S = d^T, + ie[Q, Tj\Afj, with e > 0, 
and S — > LTiR^ under chiral transformations, 

^ = '' ' ^ = ( ^ _,o/^ ) ■ (3) 

We are restricting ourselves to the two-fiavor case, but the arguments are general. In order 
for the pion field in eq. Q to be a fiuctuation about the true strong interaction ground 
state, the phases (pj are constrained so that in the expansion of eq. Q, terms linear in the 
pion field are absent. The two constraints on the phases lead to the well known relations 

rrid , iTLu 



rriu + md ' rUu + rud 



where we have used the fact that <^ 1. 



II. STRONG CP- VIOLATION IN CHIRAL PERTURBATION THEORY 

At leading order in the heavy baryon expansion, the nucleon dynamics are described by a 
Lagrange density of the form 

C = N iv ■ DN + 2gAN S^'A^ N , (5) 

where is the nucleon four-velocity and is the covariant spin operator. The chiral 
covariant derivative is given in terms of the meson vector field D^N = d^N + V^N ^ where 
= i (^t^^^ + ^^^^t). The field ^ is related to the E-field in eq.© by S = and 
^ ^ LC,W = U^R^ under chiral transformations. The leading order interaction between 
nucleons and the pions is characterized by the axial coupling constant qa ~ 1-26 in eq. (jSJ, 
where A^ = ^ (^^d^^^ — ^'''df^^^. The light quark masses contribute to the dynamics of 
nucleons through the Lagrange density 

Cm = -2aN m,5+ - 2 a iV Tr ( 771,5+ ) > (6) 

where 777,^+ = | + ^rn]^^ ^, 777, LniqR^ and 777,^+ — » Uniq^+W . The quantities a 

and a are constants that must be determined experimentally. Upon removing the GG term 
by a chiral transformation, the mass matrix becomes 777, = diag (jriue^^^ , 777^6*'^'^^ , where 0„ ^ 
are given in eq. (@]). Neglecting electromagnetic contributions to the nucleon mass splitting, 
and using light quark masses of 777^ = 5 MeV and 777^^ = 10 MeV, we find that ^ 

2a = ~ 0.26 . (9) 

777„ - 777^ 

^ The standard analysis usually invokes the approximate flavor SU{?>) symmetry, e.g. Ref. The light 
quark contribution to the baryon masses arises from 

Cm = -60 Tr [ mq^+ ]Tv[B B] - h Ti [B mg^+ B] - 62 Tr [ B B mg^+ ] , (7) 
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The sigma term is defined to be 



crjv = 2^m,— — , (10) 

where Mm is the nucleon mass in the isospin hmit m„, — >■ m, and is related to the 
quantities a and a by 

= 2(a + 2(T)m. (11) 

The value of is somewhat uncertain, with values ranging between 45 ± 8 MeV 0] and 
64 ± 7 MeV . Partially-quenched lattice computations are currently underway to evaluate 
both a and uat. 

Expanding out the interaction in eq. © to linear order in the pion field gives rise to the 
CP-violating, momentum independent NNn vertex 

£ = -^^^^iV f^°/^ 1]^]n + .... (12) 
/ mu + TUd \ n -7r''/V2 J 

It is well known that a single insertion of this interaction into a one-loop diagram gives rise 
to an electric dipole moment (edm) of the nucleon ^■ 




FIG. 1: The one-loop diagrams that contribute to the neutron edm in chiral perturbation theory. 
In QCD only vr 's participate in the loop diagram, while in partially-quenched QCD there are con- 
tributions from the bosonic mesons, M , and the fermionic mesons, x- The crossed circle denotes 
an insertion of the CP-violating vertex in eq. the square denotes an insertion of the strong 

ttNN or TT'jNN interaction from eq. |3p with derivatives promoted to covariant derivatives, and the 
small circle denotes an electromagnetic interaction with the meson from eq. ^ . 



from which one finds that, neglecting electromagnetic corrections, 

6i = — 1.1 , 52 = 2.3 , (8) 

TOs - mu ms- md 

where we have used TOj ^ 120 MeV. 
^ This set of diagrams also dominates the nucleon edm form- factor, as recently computed in Ref. [ll|. 
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The electric dipole moment of the neutron, dn, is defined by the Lagrange density de- 
scribing the interaction between a neutron and an external electric field, 

C = dnn a -Fi n , (13) 

where cr are the Pauli spin matrices, and E is an external electric field. A calculation of the 
one-loop diagrams shown in Fig. ^ leads to 

gA a e9 my, irid . ( ml\ , ^ '^^ '^d e 
2n^f^ mu + rrid \ 1^ / rUu + rud 

We have only kept the logarithmic contribution from the loop diagram, which depends 
upon the renormalization scale fi. This scale dependence is exactly compensated by the 
contributions from local counterterms, which we have combined into c(/x). There are ten 
local counterterms that contribute to the nucleon edm, as presented in Ref . , and setting 
/X ~ we anticipate that c(A^) ~ 1, where A^ ~ 1 GeV is the scale of chiral symmetry 
breaking. 

Numerically, using the value of a in eq. 0, we find the one-loop contribution to be 

dn 1.2 X 10"^*^ 9 e cm , (15) 

which is consistent with previous estimates ^ [il of the one- loop diagram . The current 
experimental upper limit is \dn\ < 6.3 x 10~^^ e cm, from which one concludes that \6\^ 5 x 
10-10. 



III. NEUTRON EDM AT FINITE VOLUME 



Lattice calculations of the neutron edm are performed on finite lattices, and therefore one 
must consider finite volume corrections in translating the lattice results to physical predic- 
tions. For lar ge e nough lattices, of course, these finite volume effects will be exponentially 
small 

0,1 

llTL 1331 1 . There has been a fair amount of work on finite volume corrections to 
quantities calculated in lattice QCD [l3,[ll,H[li[liH|2i|2l|2l|2l|2lS 
iil . 3^ 31, ll^, H^, 113, H^, but it is only recently that the properties of baryons, and in 



3il3i,E 



particular the nucleon, have been considered j3c 

In the calculations that follow, we will assume that the time direction of the lattice is 
infinite. Clearly, this can only be an approximation, but in most simulations, the time 
direction is considerably larger than the spatial directions, usually by more than a factor of 
two. By assuming that it is infinitely long, we are able to analytically perform the integral 
over energy in the loop diagrams that contribute to the observable of interest, leaving sums 
over the allowed three-momentum modes on the lattice. Details of this procedure can be 
found in Refs. [3^ , and we will not elaborate further here. By computing the one- loop 
diagrams in Fig. Q in a finite spatial volume for which the spatial dimension, L, is much 



In Refs. 



the electronic charge e is negative. 
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greater than the pion Compton wavelength, mT^L ^ 1, and for which the power counting 
rules are those of the p-regime at infinite volume, we find that 

4^) = _ Y. AV™,i|n|) , (16) 



7r2/2 niu + rud 



n^O 



where is the neutron edm at finite volume, and is its value at infinite volume. 
Kq{x) is a modified Bessel function. In Fig.|21we show the ratio d^^yd^!^^ for c(/i) = 0, as an 




2.5 3 3.5 4 4.5 5 
L(fm) 



FIG. 2: The ratio, 5 = d^ydn^\ of the neutron edm at finite volume to its value at infinite volume 
as a function of spatial lattice size L, for c{fi) = 0. 



example. The finite volume corrections are found to be quite large, primarily due to the fact 
that the leading order contribution to the edm is at the one-loop level, and not from a lower 
dimension operator. In the case of the nucleon properties previously considered js^ Is^ ]. 
such as the nucleon mass, magnetic moment and matrix elements of the axial current, the 
loop contributions are subleading, and hence the finite volume corrections are subleading in 
the effective field theory. 

As one moves into smaller volumes, where m-^j-LS, 1, the p-regime power counting is no 
longer applicable, and we move into the e'-regime 37|. In this regime, the spatial zero-modes 
are enhanced relative to the non-zero-modes, and a power counting in terms of the small 
parameter e' = m^rL is appropriate. For the neutron edm calculation, the same one-loop 
diagram shown in Fig. ^ makes the leading contribution in e' = rriTrL, and the neutron edm 
is found to be 



-2 gA a e 9 rriu mg 



+ 



(17) 



where the ellipses denote terms higher order in the e'-expansion. The classic exponential 
behavior of the p-regime becomes power law, behavior as the spatial volume decreases. 
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IV. NEUTRON EDM IN PARTIALLY-QUENCHED QCD 



While partially-quenched QCD (PQQCD) [23, 123, lid, |4l|, 142, liS, |4j is not a theory that 
describes nature, it is a theory that can be used to describe unphysical lattice calculations, 
and allows the direct extraction of QCD observables via an extrapolation in quark-masses. 
In calculating quantities in lattice QCD, the quark masses used in the generation of gauge 
field configurations does not have to be the same as the quark masses of the propagators 
computed on those configurations. The reason why this is a useful concept is that the 
computer time required to generate a dynamical configuration grows rapidly as the quark 
mass is reduced, while the time to compute a propagator grows more slowly. Currently, 
lattice calculations cannot be performed at the physical quark masses, but we wish to be 
as "close as possible" to the physical values in order to minimize the impact of quark mass 
extrapolations. 

The Lagrange density describing the quark sector of PQQCD is 

£ = [^JP~^q]1 Qn - -G^'^^'-G^^J + , (18) 

k,n=u,d,u,d,j,l 

where the left- and right-handed valence, sea, and ghost quarks are combined into column 
vectors 

Ql = (u,d,j,l,u,dj^ , Qr = (u,d,j,l,u,dj^ . (19) 

The objects r/^ correspond to the parity of the component of Qk, with rjk = +1 for k = 
1,2,3,4 and r]k = for k = 5,6. The Ql,r in eq. (|19|) transform in the fundamental 
representation of SU{4:\2)l,r respectively. The ground floor of Ql transforms as a (4, 1) of 
SU{4:)qL ® SU{2)qL whilc the first floor transforms as (1, 2), and the right handed field Qr 
transforms analogously. In the absence of quark masses, mq = 0, the Lagrange density in 
eq. (fTHj) has a graded symmetry U{4:\2)l ® U{4:\2)r, where the left- and right-handed quark 
fields transform as Q^, — > UlQl and Qr UrQr respectively. The strong anomaly reduces 
the symmetry of the theory, which can be taken to be SU{A\2)l ® SU{A\2)r ® u\l)v 0. 
It is assumed that this symmetry is spontaneously broken SU (4|2) ^ ® SU (4|2)/j ®U{l)y — *• 
5'f/(4|2)y ® U{l)v so that an identification with QCD can be made. The mass matrix, 
mg, has entries mq = disig{rnu,mii,mj,'mi,mu,md), (i.e. the valence quarks and ghosts are 
degenerate) so that the contribution to the determinant in the path integral from the valence 
quarks and ghosts exactly cancel, leaving the contribution from the sea quarks alone. This 
makes clear why the partially-quenched theory describes lattice calculations with sea quarks 
and valence quarks of differing mass. The details concerning the construction of partially- 



jerturbation theory (FQ^FT) are well known and can be found in several 
45I 46^ . The quantity that has "physical" impact for lattice calculations of 
strong CF- violating quantities is ^ 6* = — arg (sdet (mq)). 



quenched chiral 
works, e.g. Ref. 



sdet 



A b\ dct (^ - BD-^C) 



C D 



det [D) 



(20) 



The strong interaction dynamics of the pseudo-Goldstone bosons are described at leading 
order (LO) in PQxPT by a Lagrange density of the form, 



str 



+ A Vstr 
4 



mgEt + m^S 1 + «$9^$o9^$o - rnl^l , (21) 



where a$ and mo are quantities that do not vanish in the chiral hmit. In order to simply 
project out the singlet of the graded group one takes the limit mo 
field is incorporated in E via 



oo 



3- The meson 



exp 



V X 



X' 
M 



(22) 



where M and M are matrices containing bosonic mesons while x ^ire matrices con- 

taining fermionic mesons, with 



M 



( Vu 


7r+ 


JO 






Vd 


J- 






J+ 


m 


y+ 




1° 







M 



Xriu 



Xl+ 
Xlo 



(23) 



where the upper 2x2 block of M is the usual triplet plus singlet of pseudo-scalar mesons while 
the remaining entries correspond to mesons formed from the sea quarks. The convention we 
use corresponds to / ~ 132 MeV, and the charge assignments have been made using an 
electromagnetic charge matrix of Q^^'^^ = |diag (2, — 1, 2, — 1, 2, — 1). For the calculations 
we will be performing, the flavor singlet pseudo-Goldstone boson does not contribute, and 
so we do not discuss it and its associated hairpin interactions. 

The free Lagrange density describing the interactions of the nucleon and its superpartners 
which are embedded in the 70 dimensional irreducible representation of S'[/(4|2) Bijk is, at 



LO in the heavy baryon expansion [43, |4a, |49|, |5 



C = HBvVB] -2aZ^^ 



BBM+) - 213^"^^ 



BM+B] -2a£^^ 



BB) str(7W+) ,(24) 



where = \ {i'^rnq^'^ + ^^o'C)- The brackets (...) denote contraction of Lorentz and 
flavor indices as defined in Ref. (43|. 

The Lagrange density describing the interactions of the 70 with the pseudo-Goldstone 
bosons at LO in the chiral expansion is |43| . 

C = 2p {bS^^BA^) + 2p {bS^'A^b) , (25) 

where is the covariant spin vector 47, 4^ 4^. Restricting ourselves to the valence sector, 
we can compare eq. ()25|) with the LO interaction Lagrange density of QCD, 



C = 2gA NS^A^N + giNS^N tr [ A^ 



(26) 
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and find that at tree level 



P 



4 1 

-^9A + -9. 



^91 - ^9A 



(27) 



The contribution to the strong anomaly from the valence quarks is exactly canceled by 
the contribution from the ghosts. Therefore, chiral transformations of the sea quarks alone 
remove the 6'-term from the Lagrange density in eq. (fTSj) . Upon a chiral transformation 
of the valence quark, sea quark and ghost fields, the quark super-mass matrix becomes 
rriQ = diag(mue*'^", m^e*'^'', mjC*'^-'', rrt/c*'^', m^e*'^", m^e*''^'') subject to the constraint that 9 = 
— J2 (— )''*+^ (pk- The vacuum stability condition for small 6 further provides the constraint 
^u(pu = ^d4>d = = ^i4>i- Therefore, we have 



6 mi 



rrij + mi 



6 mj 



mj + mi 



9 mj mi 1 



9 mj mi 1 



mj + mi mu 



mj + mi md 



— ■ (28) 



By using this phase-rotated mass matrix in the Lagrange density of eq. (j24p . one induces 
a CP-violating interaction between the pseudo-Goldstone bosons and the baryons of the 
partially-quenched theory, in precisely the same way as in QCD. Further, this vertex gener- 
ates the leading contribution to the neutron edm through the one-loop diagrams analogous 
to those in Fig. ^ One further slight complication that can be considered is that the elec- 
tric charge matrix in the partially-quenched theory is not specified by nature; all that is 
required is that one reproduces QCD in the limit that the sea and valence quarks become 



degenerate ^ j45|, |4 
form ^ (^jg^g 



5lL l52| . In our computations, we use an electric charge matrix of the 



Working in the isospin limit where m^ 



mi 



and defining qji = qj + qi, we find 



that the leading order contribution to the neutron edm is 



e 9 msea 

e 
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A2 

X 



lo£ 



P' 



+ Fj loe 



m 



P' 



m. 



C{fi) + d (msea - ^val) + f qjl ( ^sea '^val j 



(29) 



where mj is the mass of the Goldstone boson composed of a sea quark and a valence quark, 
and 



9A 



^ (PQ) ^(PQ) ^ 



3 

1 (I3I 



(In 



Fj - 9AaM (3+yj - 9i\^ 



(30) 



^ Even this constraint is excessive. It is sufficient to determine matrix elements of operators transforming 
in the singlet and adjoint representations of the graded group. 
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As we can make the tree level identification a = {2a^^j^^ — (3^^j^^)/3, the expression in 
eq. reduces to the QCD result in eq. (fT^ when rrisea. — * ^valence and mj — > m^r, 
since F,^ + Fj = Qac^- It is important to notice that the counterterm that contributes in 
the partially-quenched case, c(/i), is the same as in the QCD case, while the other two 
counterterms, d and /, make a vanishing contribution in QCD. The expression in eq. ()29|) 
exhibits one of the well known pathologies of the partially quenched theory. One sees that 
this expression behaves as ~ msea log (mvaience)- For a fixed sea quark mass, the one-loop 
contribution diverges as the valence quarks move toward the chiral limit, in contrast to the 
case of QCD where the diagram diminishes as ~ log {ml). 

The finite volume corrections resulting from a partially-quenched calculation are obviously 
more complicated than in QCD. In the limit where the volume is large compared to the 
Compton wavelength of both the valence and sea mesons, one can use the power counting 
of the p-regime to find that 

diPQ)iL) ^ ^(PQ)(oo) _ e9m^ [FjSj + F^ 5.] 

= J2 Ko{m^L\n\) , Sj = J2 Ko{mjL\n\) . (31) 

117^0 n^O 

One can imagine performing a calculation of the neutron edm for lattice parameters such 
that myr-L <^ 1 but mjL> 1. Parametrically, we could arrange for mj^/A^ ~ e'^, A^L ~ 1/e' 
and mj/A^ e' . In such a scenario, the finite volume correction would become 

'*r'"'" = -^5^f. +■■• ■ (32) 

This somewhat bizarre computational set-up allows one to quite dramatically separate the 
contributions to the neutron edm, as the leading contribution results from one-loop graphs 
involving pions, and the contribution from mesons involving the sea quarks is suppressed. 
However, the sea quarks play a central role via the CP-violating pion-nucleon coupling. In 
the more symmetric scenario in which mT^L,mjL^ 1, the finite volume expression becomes 



p 13 



5l ^ Ll 

■:> I o 

mi rrij 

TT J _ 



+ ■■■ . (33) 



V. CONCLUSIONS 

A non-zero electric dipole moment of the neutron would provide direct evidence for time- 
reversal violation in nature. It continues to be the focus of ever more precise experimental 
measurements, and the fact that it has not been observed at the present limits of experimen- 
tal resolution provides one of the more intriguing puzzles in modern physics. In this work we 
have considered how lattice QCD calculations of the neutron edm originating from the QCD 
6'-term, performed in a finite volume and at unphysical quark masses, are related to its value 
in nature. We have provided explicit formulas that allow for the extrapolation from finite 
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volume calculations to the infinite volume limit and for the chiral extrapolation of partially- 
quenched calculations. In order for these formulas to be useful, lattice QCD calculations 
of both the light quark mass dependence of the nucleon mass, a, and the neutron edm are 
required. With the lattice value of a known with a given precision, the lattice determination 
of the neutron edm will then allow for the counterterm c(/i) to be determined. Once these 
constants are computed, the chiral extrapolation of the neutron edm to the physical quark 
masses, and to infinite volume, is possible. 
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